Behavior of magnetic impurities in gapless Fermi systems 
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In a number of systems, including certain semiconductors and unconventional superconductors, the 
effective density of states varies near the Fermi energy like \E—Ef\ t - The behavior of dilute magnetic 
impurities in such systems is studied using a non-perturbative renormalization-group approach. 
Close to particle-hole symmetry, the Kondo effect is suppressed for the cases of greatest relevance 
(r = 1 and 2). Away from this symmetry, any quenching of the impurity moment is accompanied 
by a low-temperature decrease in the impurity resistivity, rather than the increase found in metals. 
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In an interesting class of "gapless" Fermi systems, the 
density of states p(e) vanishes right at the Fermi energy 
Ep, but not at any other energy near e = E — Ep = 0. 
For instance, the quasiparticle density of states in an un- 
conventional superconductor can vary like |e| or |e| 2 near 
line or point nodes in the gap fl. Heavy- fermion and 
cuprate superconductors are candidates for this behavior. 
The valence and conduction bands of certain semiconduc- 
tors — including Pbi-^Sn^Te at a critical composition 
0, and PbTe-SnTe heterojunctions || — touch in such a 
way that, for small |e|, p(e) is proportional to |e| d_1 in d 
spatial dimensions. Electrons in a strong magnetic field 
(at least in the absence of disorder) 0] and exotic phases 
of the Hubbard model Q are also predicted to exhibit a 
linear pseudo-gap in two dimensions. 

Antiferromagnetic coupling between magnetic impu- 
rities and a metallic conduction band leads to a low- 
temperature increase in the resistivity and a reduction 
in the Curie term in the susceptibility. This Kondo ef- 
fect depends on the existence of electronic excitations 
down to zero energy, and thus cannot be fully realized in 
systems with a finite energy gap. Gapless systems with 
a density of states varying like |e| r constitute a marginal 
case, first studied in Ref. ||. Poor-man's scaling for the 
spin-i (impurity degeneracy A^ = 2) Kondo model, and 
a large- N treatment restricted to r < i (but recently ex- 
tended to include r — 1), both showed that a Kondo 
effect takes place only if the electron-impurity exchange 
J exceeds a critical value, J c cx r; otherwise, the impurity 
decouples from the band. A large- A?" treatment of mag- 
netic impurities in gapless superconductors M yielded 
similar results, except that for r < 1 or N = 2, any fi- 
nite impurity concentration was found to drive J c to zero. 
Recently, however, third-order scaling was applied to the 
n-channel Kondo model to show that, at least for n 3> 1, 
no Kondo effect can occur if r > l/(2n) 0. 

This paper reports the results of a nonperturbative 
renormalization-group (RG) study of a spin-^ impurity 
interacting with an electronic band in which p(e) takes 
one of several functional forms, each varying like |e| r near 
e = 0. A stability analysis of the RG fixed points and 
numerical calculations of impurity thermodynamic prop- 



erties are presented. Particle-hole asymmetry is identi- 
fied as a key factor in determining the low-temperature 
physics. At small asymmetry, the critical coupling J c 
above which the impurity moment is screened becomes 
so large for all r > g that the Kondo effect is suppressed. 
Larger asymmetries produce a low-temperature decrease 
in the impurity resistivity, rather than the monotonic in- 
crease which characterizes the metallic Kondo effect. 

After the completion of this work, the author be- 
came aware of a similar study, limited to particlc-holc- 
symmetric systems and a pure power-law density of states 
with < r < 1 pi. The present paper confirms the con- 
clusions of Ref. |1C| for this special case, but shows that 
the generic behavior of gapless systems is quite different. 

The Kondo Hamiltonian describing impurity potential 
and exchange scattering of band electrons can be written 

H = ffbaad + #imp = D ^ dee c\ a c erT 

a 

+Vj2flf0* + ■'E/.U'W/O.r' ' S - (1) 
cr cr, a' 

Here e = (E — Ep)/D is a reduced kinetic energy, mea- 
sured from the Fermi level in an isotropic band of width 
2D; J > represents antiferromagnetic exchange; c ea an- 
nihilates a spin-cr electron in an s-wave state of reduced 
energy e; and f 0a = Jde^/ Dp(eD) c ea is the combina- 
tion of c £(T 's that destroys an electron at the impurity 
site. The fcrmionic operators are normalized such that 



{4a: c e , ff ,} = S{e - e')<W and f 0a ,} 



Gapless systems can be modeled by a density of states 



p{e) = p(eD) 



Po\e\ 




if |e| < 1, 
otherwise, 



(2) 



with r > 0. This form is oversimplified, but it will 
turn out that for practical purposes the results arc lit- 
tle changed if one introduces band asymmetry, limits the 
power-law variation to the vicinity of Ep, or allows p(0) 
to be small but nonzero. Standard results |L2j should 
(and will) be recovered in the metallic limit, r = 0. 

Equation ([l]) can be treated nonperturbatively using 
a generalization of Wilson's numerical RG method [12] 
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to an arbitrary density of states. The Hamiltonian is 
written H = liniM-^oo DK~ m / 2 Hm, where 

H M = A 1/2 H M -1 + ]T [smfLfMa 



+ 6f-i (/m ct /m-i, ct + h.c.) 



(3) 



for M > 0. Here H-i = H imp /(A^ 2 D) and £_ x = 0; 
A > 1 parameterizes the logarithmic discretization of the 
conduction band Jl2[ ; and /m annihilates an electron in 
a state, centered on the impurity site, that becomes pro- 
gressively more delocalized as M increases. Numerical 
solutions of successive Hm's capture the key physics at 
a sequence of temperatures T « DA~ M / 2 p^Jl^ |. 

All dependence on p(e) is contained in £m and em, 
which are given by a set of Lanczos recursion relations. 
Quite generally, if p(e) = p(—s), then em = for all M. 
For the density of states in Eq. (||) , one finds [|l4| that 



M- 



1+r l-A-< 2+r ) 



M 



2+r l-A-(!+ r ) 



A~ 1/2 A~ (JU mod 2)r/2 . (4) 



Preliminary analysis: Even before one resorts to heavy 
computation, study of the limits J — and oo provides 
valuable insight concerning the density of states given by 
Eq. (||): (1) For all r > 0, there is a weak-coupling regime 
in which the Kondo interaction is irrelevant ||. (2) At 
particle- hole symmetry, a Kondo effect (i.e., scaling to 
infinite J) is ruled out for all r > | (3) If V + 

in Eq. (|l|) , there is a stable strong-coupling limit in which 
the impurity is completely screened, but the impurity re- 
sistivity is zero rather than taking its maximum possible 
value, as it does in metals. (4) Since each RG trajectory 
must flow from an unstable fixed point to a stable one, the 
intermediate-coupling region is also constrained. If the 
limit J = oo is stable, at least one (unstable) fixed point 
must lie at some < J c < oo, whereas instability about 
J = oo is consistent with uninterrupted RG flow from 
strong to weak coupling. The derivation of these results, 
outlined in the next four paragraphs, follows Ref. 113]. 
For J = 0, different electron spins decouple. Each spin 

AO) 



Simp = In 2 p3| . For r > 0, electrons at the Fermi energy 
scatter from the impurity with a phase shift 5(e = 0) = 0. 

Any deviation of H M from H M \v) must be a com- 
bination of /t's and /'s which respects all symmetries 
of Eqs. (Q) and (||), multiplied by the same factor of 
A M / 2 as enters Eq. (||). Using the above expansion of 
the /'s and the relation T « DA~ M / 2 , one can clas- 
sify the relevance of any such operator |l2|. For in- 
stance, the exchange and potential-scattering operators, 
Oj = A M l 2 f\ a \a- aa i f Qa , ■ S and O v = A M / 2 flf 0tT , both 
vary like T r , a temperature dependence which can be ab- 
sorbed into effective couplings J e e(T) = J • (T/D) r and 
V e s(T) — V ■ (T/D) r . Any other allowed operator has an 
algebraically greater temperature exponent. Thus, for 
all r > 0, Oj and Ov are the leading irrelevant opera- 
tors about a stable J = fixed point. In the metallic 
case (r = 0), by contrast, J c g grows logarithmically as T 
decreases, so the fixed point is marginally unstable. 

In the limit J = oo (with V finite), an f electron is 
locked into a spin singlet with the impurity. Hopping to 
or from / states is ruled out, so different electron spins 
again decouple. The effective Hamiltonian is H^(0) 
given by Eq. M). Now f m oc A-( 1 ~ ri ) M / 4 for m odd and 
f m oc A~ ( - 3 ~ r ^ M / 4 for m even, where r n = min(r, n). 
One finds that T^imp = r i/8, S^mp = 2riln2, and 
<5(0) = (I — ri)n/2, indicating that even at J = oo the 
impurity degree of freedom is not completely quenched, 
but instead is partially transferred to the band [0 . 

The stability of this fixed point hinges on whether or 
not the problem is particlc-holc symmetric. If V = 
in Eq. (]]]), the most relevant perturbation about Hj^ 
is A M/2 {flf 1 - \) 2 oc T 1 - 2 '' 1 , in which case the fixed 
point is stable for r < |, but is destabilized for all 
r > h- Setting V admits an additional operator, 

A M/2 /L/i<r oc T- ri . For r > 0, this freezes the /i states 
as T — > and drives the system to a new J — oo fixed 
point, described by H M (0). Due to the "odd-even" char- 



acter of the Hamiltonians H 



(L) M 2 ) 



M J 



H 



(L) 
M 



H\ 4 '(0) has the same 
low-energy spectrum and stability properties as the J = 
Hamiltonian H^(0). At this fixed point, T^imp = 0, 
is described by an effective Hamiltonian H$(V), where S imp = 0, and 6(0) = ir; the impurity contribution to the 

resistivity, g oc sin 2 6(0), vanishes for r > 0. 

Numerical results: The RG picture has been completed 
by solving the Hamiltonians Hm in Eq. (||) for arbitrary 
J. The entropy, heat capacity and magnetic susceptibil- 
ity have been computed by adapting methods developed 
for the metallic Kondo problem p2lff5|. (Space permits 
only the susceptibility to be plotted.) The main sources 
of error are the discretization of the conduction band and 
the truncation of the basis of Hm- In this work 800-2000 
states were retained, sufficiently many that truncation er- 
rors are negligible (smaller than the symbols in Figs. |l|- 
||). All data shown were obtained using a discretization 
parameter A = 3. Changing A affects the results only 



(V)=A M ' 2 [(V/D)flf Q 
M-l 

+ J2 A-" l / 2 £ m (/tj m+1 + h.c.) 



711 — L 



(5) 



For M S> 1, it is found numerically that the low- lying 
eigenvalues of (V) are independent of M and V, 
while the /'s have simple expansions in terms of the ex- 
act cigenoperators go,..., 9m- fo = A~( 1+r ) M / 4 £\ a^gj, 
where ccj is independent of M; all other / m 's decay with 
increasing M at least as fast as fo decays. The impurity 
susceptibility is Curie-like, T^imp = j, and the entropy is 
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FIG. 1. Critical Kondo coupling J c /D vs. the power r 
entering Eq. (^), for potential scatterings V/D = 0, 0.1 and 
0.5. (The connecting lines are provided as a guide to the eye.) 
For r = 1 only, data are shown for additional values of V. 

through small shifts in the value of J c . 

Results are presented first for pure power-law densities 
of states [Eq. (ph], confirming and extending the conclu- 
sions drawn above. Consider the particle-hole symmetric 
case, V = 0. For < r < i there exists a critical cou- 
pling J c (r), plotted in Fig. pi. For small r, J c w 2rD (in 
agreement with Ref. 0), but the curve J c {r) turns up- 
ward and then terminates at r = | . Far below a crossover 
temperature, Tx ~ — J c \/ max( J, J c )] 1/,r , any cou- 
pling J <J C yields the same excitation spectrum as , 
implying renormalization to zero exchange; whereas val- 
ues J > J c reproduce the strong-coupling spectrum of 
Hj^. In the region r > 5, the critical point disappears, 
so any exchange < J < 00 renormalizes to zero . 

Turning to the case V ^ 0, one finds as expected that 
the curve J c (r,V) extends to arbitrary r (see Fig. [l]). 
For r > i, the critical coupling is strongly V^-dcpcndcnt: 
As \V\ increases from zero, J c initially drops from in- 
finity. However, once the potential scattering becomes 
large enough to inhibit hopping of electrons to or from 
fo states, further increases in \V\ serve only to push J c 
higher. This is illustrated by the r = 1 data in Fig. [l]. 

Figure ^| compares impurity susceptibilities for r = 0, 
0.3 and 0.6. For r = 0, Tximp falls monotonically to zero 
as the temperature decreases — a sign of the complete 
screening of the impurity moment which occurs for any 
J > 0. This should be contrasted with the r = 0.3 data: 
For J — 0.3-D, T^jmp rises at low temperatures towards 
the weak-coupling value 4; whereas the J = 0.57D curve 
is almost flat, indicating that J is very close to J c . Nei- 
ther curve is greatly affected by potential scattering (not 
shown in Fig. ||), and in both cases Simp ~ In 2 over the 
entire temperature range shown. The two J = 0.8D 
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FIG. 2. Impurity susceptibility T^imp vs. logfT/D), for 
three power-law densities of states denned in Eq. (|2|). 

curves do exhibit impurity screening at temperatures 
T < Tx ~ 0.05D. For V = 0, the susceptibility remains 
Curie-like, reaching a limit Tximp = 0.0377 very close to 
the strong-coupling value r/8. The computed entropy, 
Simp = 0.416, is also in good agreement with the pre- 
dicted value 2rln2. Setting V — 0.1D causes a second 
crossover, around T' x = \V/D\ x l r T x ~ 2xlO" 5 L», to the 
particle- hole- asymmetric fixed point, at which T^imp = 
and S'imp = 0. Note, though, that the V = fixed point 
still dominates the behavior in the range T' x < T < Tx- 

Figure || also shows two r = 0.6 curves. For V = 0, 
the impurity remains unscreened at low temperatures, 
as expected from the vanishing of the intermediate fixed 
point. If instead V = 0.1D, there is a finite critical cou- 
pling, J c w 1.3-D; now an exchange J = l.bD ensures 
that the impurity is screened for T <C T' x « 10~ 3 D. 

In order to investigate the sensitivity of the preceding 
results to the precise form of Eq. (^) , the effect of three 
changes to the density of states will be briefly examined: 

(i) In real systems, the power-law variation of p(e) is 
restricted to an energy range |e| < A, with p(e) ss p(A) 
for A < |e| < D. At temperatures T ^> A, one finds the 
standard Kondo physics of metallic systems. For T <C A, 
however, the impurity "sees" a pure power-law density of 
states, so the results above still apply, albeit with J re- 
placed by J ff(A) > J. This enhancement makes realiza- 
tion of the Kondo effect more plausible for systems with 
r > 1, a range in which the critical couplings would other- 
wise be unphysically large. Figure [| illustrates the point: 
with A = 10" 3 D and V = 0.3D, a coupling J = 0.5D 
fully screens the impurity, for both r = 1 and 2. 

(ii) Band asymmetry can be introduced by shifting the 
Fermi level away from the band center, while retaining 
the power-law variation of p(e) about the new Fermi en- 
ergy. Although this changes the £m— l's entering Eq. (||), 
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T/D 

FIG. 3. Impurity susceptibility, Tximp vs. log(T/D), 
for restricted power-law densities of states in which 
p(e) =po|A/D| r for A < |e| < D. The r = oo curves rep- 
resent an insulator. 



a sufficiently large exchange coupling, J > J c , results in 
complete screening of the impurity, marked by the van- 
ishing of Tximp and a low-temperature decrease in the 
impurity resistivity; for J < J c , the impurity asymp- 
totically decouples from the band. Systems which lie 
close to particle-hole symmetry are unlikely to exhibit a 
Kondo effect because in such cases J c becomes unphysi- 
cally large. Further work is needed to extend these results 
to degeneracies N > 2, to obtain the full temperature de- 
pendence of the resistivity, and to treat self-consistcntly 
magnetic impurities in unconventional superconductors. 

I thank P. Hirschfeld and A. Schiller for many helpful 
discussions; C. Cassanello, E. Fradkin, and K. Ziegler for 
informative conversations; L. Oliveira for pointing out 
the form of Eq. (0); and T. Costi for drawing my atten- 
tion to Rcf. |lQ]. This work was supported in part by 
NSF Grant No. DMR-9316587. 



and generates nonzero coefficients e m , the physical effects 
are little different from those of potential scattering. 

(iii) Partial filling of the pseudo-gap may introduce a 
lower cutoff A' on the power-law density of states, with 
p(e) w p(A') for |e| < A'. The finite value of p(0) must 
eventually produce a standard, metallic Kondo effect. 
However, for r > 1 and most plausible values of J, A 
and A', this occurs at inaccessibly low temperatures. 

Discussion: Although large- A studies of gapless sys- 
tems have treated symmetric bands with zero po- 
tential scattering, they have found no sign of the Kondo 
effect disappearing for any power r. This discrepancy 
with Ref. jlC| ] and the present work may stem from the 
mean-field nature of the large-TV method, or from the 
symmetry-breaking that is implicit, for all N > 2, in the 
restriction that the impurity level be singly occupied. 

The r = oo limit of the restricted density of states 
defined in (i) above describes an insulator with gap 2A. 
The results of this work are entirely consistent with those 
known for gapped systems. At particle-hole symmetry, 
an impurity in an insulator retains its moment, no matter 
how large J is made; away from this symmetry, the spin 
is screened provided that J > J c « 2D/ln(L»/A) [|§. 
The mapping of an impurity in an s-wave BCS supercon- 
ductor onto Eq. (^) necessarily introduces particle-hole 
asymmetry, resulting in a finite J c p7|| . Figure ^ demon- 
strates the similarities between the impurity susceptibil- 
ities for r = 1, 2 and oo. However, the low-temperature 
decrease in the impurity resistivity is a signature of gap- 
less systems which has no counterpart in insulators. 

This study shows that the behavior of dilute spin impu- 
rities in gapless Fermi systems differs qualitatively from 
that in metals or insulators. In most gapless materials, 
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